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On  Sums  of  Lognormal  Random  Variables* 
by 
E.  Barouch  and  Gordon  M.  Kaufman 


1 .   Introduction 

The  lognormal  distribution  has  been  used  as  a  model  for  empirical 
data  generating  processes  in  a  wide  variety  of  disciplines.   Aitcheson 
and  Brown  [  1  ]  cite  over  100  applications.   In  portfolio  analysis 
(Lintner  [  3 ])  and  in  statistical  studies  of  the  deposition  of  mineral 
resources  (Barouch  and  Kaufman  [  2  ]  and  Uhler  [  4  ])  sums  X  +...+X   =  K 
of  lognormal  random  variables  (rvs)  are  of  central  interest. 

Here  we  compute  approximations  to  the  density  of  a  sum  of  N 
mutually  independent  and  identically  distributed  lognoirmal  rvs.   As  is 
well  known,  the  density  of  a  sum  of  independent,  identically  distributed 
rvs  is  given  by  the  inverse  Fourier  (or  LaPlace)  transform  of  the  Nth 
power  of  the  characteristic  function,  so  we  begin  by  studying  the 

characteristic  function  of  a  lognormal  density.   We  perform  an  asymptotic 

2 
analysis  of  it  in  its  various  regions  for  N  and  Var(X.)  =  a      large  and 

compute  approximations  to  the  density  of  the  sum  by  transforming  back. 

We  find  that  (a)  for  values  of  K  larger  than  its  mean,  the 

density  of  K  is  approximately  a  three  parameter  lognormal  density 

(cf.  (44));  (b)  for  values  of  K  near  its  mean,  the  densitv  contains 

both  lognormal-like  and  normal-like  components,  (cf.(45))-  (c)  for 

values  of  K  larger  than  order  one  but  smaller  than  its  mean,  the  density 

is  approximately  a  three  parameter  lognormal  density,  (cf.(46))  and  (d)  for 

values  of  K  smaller  than  order  one,  the  density  is  approximately  lognormal 

(cf.(47)).  072G74I   ■ 

Tf~TK^authors  thank  their  colleague  Hung  Cheng  for  a  very  fruitful  discussion. 


2.   Lognormal  Characteristic  Function 

2. 1  Properties  of  the  Characteristic  Function  and  Approximations 

The  characteristic  function  G(y)  of  the  lognormal  distribution 
is  defined  by 

G(y)  =  /  exp{-iyx ^^^S  x}-^  (1) 

where  y  is  complex  with  Imy  <_   0  and  y   has  been  set  equal  to  0  without 
loss  of  generality.   The  function  G(y)  is  analytic  everywhere  in  the 
lower  half  of  the  complex  y  plane,  and  continuous  from  below  for  real  y. 
However,  G(y)  is  not  analytic  near  y  =  0,  and  this  greatly  enhances  the 

difficulty  of  approximating  G(y). 

""ivx 
An  obvious  way  to  attempt  computation  of  G(y)  is  to  expand  e 

in  a  Taylor  series  around  zero  and  integrate  term  by  term;  i.e.  to  . 

express  G(y)  in  a  moment  series.   In  so  doing,  we  are  expanding  G(y) 

around  a  point  at  which  G(y)  is  non-analytic.   Hence  it  is  not  at.  all 

surprising  that  this  expansion  fails  in  every  respect.   The  first  few 

terms  of  such  an  expansion  cannot  be  looked  upon  as  a  "small  y"  approximation 

to  G(y)  as  the  resulting  polynomial  is  analytic  while  G(y)  is  not. 

1  2  2,„ 
T,,  .  T/.n\  .  ■  .  ./.\nnian/<i 
Furthermore,  since  the  (n+l)st  term  m  the  series  is  (-i)  y  — j-e      , 

this  moment  series  is  divergent  for  all  y  ^  0.   Trying  hard,  the  moment 

2 
series  can  be  viewed  as  an  asymptotic  series  provided  that  o   is  very 

2 

small.   This  fact  is  not  very  useful,  since  for  a     small  enough,  the 

—iy 
series  is  well  approximated  by  e    ,  the  characteristic  function  of  a 

density  concentrated  on  the  point  1. 


Consider  a  point  y  =  5-in,ri>0  and  rewrite  G(y)  as 

oo 

G(y)  =  — —     /  exp{-nx  -  i^x  -  -^log^x}^         (2) 
a/li^       0  2a       ^ 

-i£x 
As  long  as  n  remains  positive  we  may  expand  e     in  a  power  series, 

since  we  are  expanding  G(y)  around  a  point  in  its  domain  of  analyticity. 

Hence  we  may  write 

00  j  oo 

G(y)  =  I     -^^TT —     /  exp{-nx  -  ^log^xlx^'^dx     (3) 

j-=0   ■^  ■    a/27   0         2a 

The  limit  n  ^  0  is  not  allowed  after  expanding  since  these  two  operations 
do  not  commute  and  a  study  of  G(y)  based  on  (3)  must  be  done  with  extreme 
caution.   Each  term  in  (3)  can  be  derived  from  the  first  term,  by 
differentiating  with  respect  to  n-   This  is  allowed  since  the  series 
is  uniformly  convergent. 

Consequently,  to  construct  a  good  approximation  to  G(y)  we  need 
only  study  its  behavior  for  y  lying  on  the  negative  imaginary  axis.   This 
of  course  is  not  surprising,  since  any  operation  on  G(y)  can  be  performed 
by  deforming  the  contour  of  integration  through  the  axis  Imy  <  0. 

We  now  focus  Our  attention  on 

00 

G(y)  =  /  exp{-yx rlog  x] 

a/27  0         2a       ^ 

for  real  positive  y  (namely  y  =  -in,  n  >  0,  and  so  henceforth  n  =  y) . 

2 
We  further  assume  that  a     is  large. 


Immediate  properties  of  G(y)  are: 

(i)  G(0)  =  1 
(ii)   lim  G(y)  =  0 

y->oo 

(iii)   |G(y)|  <_   1 

2       2 

(iv)  -—  =  -e    G(ye  ) 

dy 

2       2      2 
,  ,   3  G    2a  „.   2a  , 

8y 

From  properties  (iv)  and  (v) ,  it  is  apparent  that  differentiating 
2        2 
G(y)  rescales  y  by  e   .   If  e   is  large,  even  if  y  is  small,  sufficient 

differentiation  moves  G(y)  to  its  asymptotic  region.   For  example  if  y  =  0(1), 

2 

-r—  is  a  constant  times  G(ye   ).   Thus,  when  we  wish  to  perform  an- integration 

°y 

involving  G(y),  asjonptotic  evaluation  of  G(y)  may  be  necessary. 

2 
Let  y  ^  °°,  and  log  y/a  be  large.   Change  variables  in  G(y):   yx  =  z 

to  obtain 

G(y)  =  exp{-  ^log^y}  /  exp{-z  -  -iylog^z}z^   1°8  y  dz 

a/lt  2a  0        2a  z    ^  ^ 

Since  — ^-r-^   is  assumed  large  and  positive,  the  major  contribution 

a  12 

to  (4)  cannot  come  from  z  ~  0.   Thus,  — r-  log  z  is  small  compared  with 

2a 

z,  and  can  be  dropped.   This  is  of  course  a  crude  approximation.   To 

1     2 
improve  it  we  expand  exp{-  — j   log  z}  in  a  power  series.   As  long  as 

2  2a' 

log  y/a  >  0  the  resulting  term  by  term  integration  is  uniformly  convergent. 

Equation  (4)  takes  the  form 


_  1  j^  2 

a/27r  j=0   2a    j!  -"q  ^         z 

(5) 


The  integrals  in  (5)  can  be  expressed  in  terms  of  derivatives  of 
the  r  function: 

/V^log  z)2j    ,^"^logydz    =  r^'-^^a-^log  y)      (6) 


0 


with  i^  >  0. 

a 


z 


Then  G(y)  takes  the  form 

G(y)  =  — i—  exp{-  ^Wy}  I    (-  -^) ^-^^ ^^^  ^  (o~^log  y)   (7) 
c.^      2a      j=0   2a   ^• 

This  expansion  is  exact  but  not  really  useful  for  computation  since 

higher  derivatives  of  the  r  function  are  very  complicated.   We  now  make 

2 
use  of  the  assumption  that  log  y/a  is  large,  and  approximate 


o        2"   -2         -2 
r   (a'^log  y)-  ^  '^(^     log  y)r(a  log  y)  (8) 

where  4^  is  the  logarithmic  derivative  of  the  r  function.   To  justify 
this  approximation  we  write  a  table: 

r(a) 

r'  (a)  =  ilia)    r(a) 

r"(a)  =  {/(a)   +   i|;'(a)}r(a) 

r'"(a)  =  {ij/^Ca)  +  3i;;(a)i|j'(a)  +  <J;"(a)}r(a) 

r""(a)  =  {/(a)  +  6i>^ia)i>'(a)   +   44^(a)i|-"(o)  +  3ii;'^(a)  +  (p'"(a)}r(a) 


2 

with  a   =  log  y/a  .   If  we  substitute  the  leading  term  in  the  asymptotic 

series  of  \J^(a)  which  is  log  a,    il/'Ca)  ~  — ,  we  can  drop  all  derivatives 

a 

of  ij)  compared  with  i(;,  and  T     "*  (a)  -  [i|;(a)]  -^rCa).   Substitution  of  (8) 
into  (7)  yields  a  summable  series  for  G(y) : 

1         12-''            12-2 
G(y)  :=:  exp{ ylog  y}r(a  'log  y)exp{ ~  tJ;  (a  log  y)  }    (9) 

o/StT       2a^  2o 

We  can  compute  corrections  to  any  order  desired.   For  instance,  in 
order  to  obtain  the  next  term  we  write 

r^^^^a)  =  /ha)    r(a)  +  j  (2j-l)i|; '  (a)^^^"^(a)r(a)         (10) 

which  yields 

1         12-2            12-2 
G(y)  i  . — exp{ 2l°8  y}'^(°     log  y)exp{ ^    (a  log  y)} 

o-ZIt;  2a  2a 

*  (1  --i-^.(i£Y)[i  -^M>\^)]  +o([r(^)]h      (11) 

2a      a       a      a  a 

This  expansion  breaks  down  when  i^' (a)  can  no  longer  be  neglected 

relative  to  ij;(a)  and  is  particularly  bad  for  a  •''  0  (y  ~  1) . 

2 

Also,  (9)  is  valid  for  y  >  e°    ;  surprisingly,  (9)  works  well  for 
y  close  to  expi-ro  }. 

We  have  already  argued  why  y  ^  0  is  not  legitimate  here.   Properties 
(ii)  and  (iii)  are  immediate  for  (9)^  and  it  satisfies  (iv)  to  the  order 
of  approximation  computed.  [If  (iv)  is  obeyed,  the  rest  of  the  derivatives 
are  immediate].   To  see  this,  we  neglect  ^' ,   differentiate  (9),  and  show 


1  2  2 

that  the  result  equals  -exp{-ra  }  times  (9)  with  argument  ye°  in  place 

of  y.   Property  (iv)  takes  the  form 


1  12-2  12-2 

exp{-  — jlog  y}r(a      log  y)exp{ ^    (a      log  y)} 


a/lTi  2a  2a' 


2-1  2-1        -2 

X    (ya    )      log  y  +   (yo    )      i|;(a      log  y) 


4  -1        -2  -2 

-   (yo   )     ^(a     log  y)4''(a     log  y)} 


(12) 


2 


=  -    -jn  exp{-  -ijlog^y)  ^^  Ho   ^log  y)exp{-  -^^[^'(o    ^og  y)    +  t-^— ]^ 
a/27  2a  yo^  2a  ^°S  y      ^ 

Approximating 

2 

r  1      r  ,   /    ~2,  .0.2, 

exp{ 2^^(o     log  y)   +  ^  ]    } 

-2 
/   1  .|2^  -2^     .,_    \p(a     log  y), 
-  expi 2^  (a  log  y)}[l  -  ^=^-=^ &-J-1-] 

2a  log  y 

and  neglecting  4*'  ii^  the  LHS  of  (12)  we  obtain  (after  some  simplification) 
the  identity 


^  -  ^  ^(a-^log  y)  =  l^Y  ^1  -  liT'^(^>>         (1^) 
a     a  a  ^    a 


Thus,  equation  (iv)  is  obeyed  by  (9)  when  we  keep  all  large  and 
order  1  terms.   The  complicated  way  in  which  the  equation  is  obeyed  is 
due"  to  the  rich  structure  of  G(y) ,  despite  its  "simple" 
integral  representation 


Next  we  study  G(y)  for  small  y.   It  seems  reasonable  to  expect 
that  for  very  small  y,  G(y)  can  be  approximated  by  a  polynomial  composed 
of  the  first  few  terms  of  its  moment  expansion.   However,  it  must  be 
done  carefully:   since  G(y)  is  not  analytic  in  a  neighborhood  of  y  =  0 
a  meaningful  small  y  approximation  must  possess  this  feature.   Furthermore, 
since  the  asymptotic  expansion  of  G(y)  for  large  y  does  not  exhibit  an 
additive  polynomial,  one  must  show  how  the  polynomial  disappears  as  y 

increases. 

3  2 
We  begin  with  y  =  0(exp{-  -ja   }) ,  since  this  is  a  relevant  order 

of  y  for  sampling  without  replacement  and  proportional  to  random  size 

(cf.  [  2])  and  then  generalize.   Adding  and  subtracting  1  -  yx  from 

3  2 
exp{-yx}  write  G(y)  for  y  =  0(exp{-  -ra  })  as 

00 

G(y)  =  /  exp{-yx 2^°8  ^'^~^ 

a/2T\       0         2a 


(15) 


1        A   2, 
=  1  -  y  exp{-20  } 


00  -2 

J    1      r  __1  T   2-,c,-z.    ,    1    2,  a   log  y  dz 
H exp{ rlog  y}  J  (e  -l+z)exp{ jlog  z}z     ^  ^   — 

a/27      2a        0  2a  ^ 


Treating  the  above  integral  in  the  same  way  as  (4),  we  obtain 

CO  -1 

C(y)  =  1  -  yexp{^  }+  — ;:;:2sxp{ ^log  y}  I   {( ^)  -rj- 

(16) 


av2Tr      2o       j=0   2a   -' ' 


X  [-4t  /"(e-^-l+z)-!^]} 
da^        0 


for  a  =  a  log  y 


In  (16),  -2  <  a  <  -1,  and  so  the  integral  exists.   In  fact,  this 
integral  is  an  integral  representation  of  r(a)  for  a  negative  and 
non-integer  valued.   Given  (16)  the  method  used  to  compute  an  expansion 
for  large  y  can  be  used  and  we  find  that  for  y  =  0(exp{-(in  +  -7)0  })  and 
integer  m,  we  can  add  and  subtract  a  polynomial  of  the  m   degree  from 
exp{-yx}  to  obtain 

G(y)  -  I  Szlll^^pih^ah   +  -^  exp{-  -i-log^ylp (l^)exp{-  -^^ (l2M^) 
j=0  J-  a/27      2a  a^  2o^  a^ 


(17) 


for   -m  >   °^^y  >  -(m+1) 
a 


This  expansion  possesses  both  of  the  features  needed  for  it  to  be  a 
meaningful  approximation  to  G(y)  for  small  y;  it  is  non- analytic  at 

y  =  0  and  has  a  polynomial  piece.   However,  it  has  two  major 

2 
defects:   it  is  not  defined  at  y  =  exp{-ma  }  for  integer  m, 

2 
and  it  appears  as  if  G(y)  is  discontinuous  at  y  =  exp{-ma  }.   Hence 

formula  (17)  can  be  regarded  at  best  as  an  approximation  of  limited 

validity.   Furthermore,  it  does  not  explain  the  relation  between  the 

rising  power  of  the  polynomial  term  as  y  decreases  and  the  lognormal 

term.   Therefore  further  analysis  is  needed. 

We  split  the  integral  representation  of  G(y)  into  a  sum  I^  +  I 

of  two  integrals  defined  by 

I  (y)  =  J    expi-yx  -  — ^log  xj-^  and 

^     a/2?r   0  2a  ^ 

(18) 

00 

T  /  ^      If      /       1  -,   2  idx 
19(7)  =  — =1  J   exp{-yx  -  — ^log  xj— - 

^      av^   1/y         2a 


10 


Wh 


en  y  is  sufficiently  small,  I^  is  small  compared  with  Land  in  the 


limit  y  ->■  0,  I^  =  G(0)  =  1.   Consequently,  we  first  analyze  I   for  y 
small.   Expanding  exp{-yx}  and  integrating  term  by  term  we  obtain  a 
uniformly  convergent  series: 

I,  (y)    =  -4=     I     -^     /^"^expf-  -i^log'xlx^   ^  (19) 

a/27  j=0       ^"0  2a 

CO  •  -2 

1      V      (-y)-*           rl-2   2,         .    flogCye^^   ), 
=  J     I  .\      exp{yj    a    }   erfc[ — ^'-^ ] 

j=0   ^-  ^o 

with 

oo 

erfc(z)  =  -^  /  exp{-t^}dt  Hi-  erf (z)  (20) 

for  z  _>  0.   In  (19)  we  must  distinguish  between  three  types  of  terms: 

2 
setting  y  =  exp{-Xa  },  with  X  >  0 

c(j  -  X)<  0,  (21a) 

a(j  -  X)>  0,  (21b) 

a(j  -  A)-  0.  (21c) 

The  argument  in  (19)  is  — a(j  -  A),  hence  there  are  a  finite  number  of 

^  2 

terms  for  which  j  -  A  <  0.   Since  a  is  assumed  large,  there  may 

be  one   term  for  which  a(j-A)-  0  (when  j  -  A  =  0(a  ),  there  is  only 
one  such  term).   Remaining  terms  are  of  type  a(j-A)>  0. 

We  replace  erfc(*)  in  each  term  for  which  j  -  A  >  0  by  its  asymp- 
totic expansion,  namely, 

err>-(z)  =  -^  exp{-z^}[l  +  y  (_i)  ^  il^ciUl]  (22) 

^"^  £=1       (2z^)^ 


12 


"   C-1)       -2      -1 
Since  J      - — r^— [j  +  a      log  y]    is  a  series  representation  of  the 

j=0   ^"         _2 
incomplete  gamma  function  y (a   log  y;  1),  (25)  may  be  rewritten  as 


Ii  (y)  -  I     ■  I   exp{-rj  a  }  + exp{-  — j^oS  yH(^  lo8  YJ  1) 

-^  j=0        J-  '^  o/ZtT  2a 

(26) 

(-l)^"^-*-          ,        1   ,      2    ,rl        /    1   T      2,      (k+l)a^-,      ^    .   1   1  (k+l)a^, 

+  -^^ exp{ rlog  y}[-exp{ — rlog    (ye^  )}erfc( ^log  ye  ) 

(k+D!        2a  2a  /2a 

-  (a/27(l  +  k  +  a"^log  y))""*-] 

-2  -2 

Even  though  yia     log  y;  1)  has  a  pole  for  a  log  y  a  negative  integer, 

2  -2 

when  y  =  exp{-(k+l)a  }  the  pole  of  y(a      log  y;  1)  is  cancelled  by  that 

-2      -1 
of  (1  +  k  +  a  log  y)    so  that  (26)  is  a  legitimate  representation  of 

I  (y)  for  all  small  y. 
^  .    2 

We  now  turn  to  y  -  1  and  first  consider  y  >  1.   Then  log (ye   )  >  0 

or  j  -  A  >  0  for  all  j  except  j  =  0,  so 


f   1  1   2  ,  ^   (-1)^     ,  1  T   2,  ja^,,   .  ,log(ye^ 
>{ p-log  y}  I  .,   exp{ — ^log  (ye-J   )}erfc( — ^-^— 

2o       3=1   ^'      2a  i/2a 


I^(y)  =  1  exp(   '   --'-      ^  ^^^  ...-.JL_.,-2.-Ja\.,„.,.log(ye^-) 


1  1  (27) 
.1  ^  /log  y,, 

+  -r€rfc( ^-^) 

2  »^a 


The  approximation  (26)  to  I  (y)  for  small  y  was  computed  keeping  only  the 

first  term  of  the  asymptotic  series  for  erfc(');  it  is  not  clear  a  priori 

that  this  leads  to  an  accurate  approximation  to  I, (y)  when  y  =  0(1),  so 

we  replace  the  asymptotic  series  for  erfc( log  y)  with  its  series 

/2a 
expansion  for  small  argument  z  >  0, 


13 


_   2  00  2j+l 

erf(z)   =   e  J     — r ,      erfc(z)    =   1   -  erf(z), 

j=0   r(|  +  j) 

and  approximate   I    (y)   with 

I.  (y) exp{ rlog  y}    )_     — 77-   [j    +  a      log  y] 

2a/2^  2a  j=2       ^• 

1                                 1            9         °°              1  T  2j+l 

+  i[l   -   exp{-  -^log^y}    I     — r^^ (12^^)  ]  (28) 

2a  2=0   r(|  +  j)   ,^0 

2 

-  -ry  exp{-;rti  }erfc( ^log  ye   ) 

2  ■ 

When  a  is  large  the  last  term  in  (28)  can  be  incorporated  into  the 

first  sum. 

2  2 

For  y  -  1,  y  <  1,  and  a     large  enough  so  that  log(ye^  )  >  0, 

00      j  2  2 

I  (y)  =  iexp{-  -^log^y  J     izil_  exp{  1  log2(ygja  )  }erfc(-^log(yeJ^  )) 
2a  j=l   -■■      2a  /2a 

+  |[2  -  erfc(-  -^  log  y)  ] 
/2a 

=  -^  exp{-  -^log^y}  f     -^=if-   [j   +  a'^log  y]"^  (29) 

2/2TTa  2a  j=2       ^" 

1                           1          o       °°            1  1  2j+l 

+  hi  -  exp{-  ^log^y     ^     — 1 (-  1°^)  ] 

2a  j=0  r(|  +  j)  /2a 

2 

1  ,.1  2,      -    .  1     .  a  . 

-  -p'  expi^tJ    }ei:fc( log  ye      ) 

/2a 
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As  y  -»■  1 ,  both  (28)  and  (29)  approach  the  same  limit,  namely. 


1    y  Izili 


lim  I,(y)  -  i{l  -  exp{4j^}erfc(-2-)  + 

y-^1  ^      '^  /I         a/27  j=2 

(30) 

=  i{l  -  exp{^^}erfc(-^)  +  —^[1  +  Ei(-l)  -  y]} 

/2    0/2^ 

where  y  =  .57721  56649,  Euler's  constant  and 

Ei(-l)  =  -  /  e'^^t"  dt  -  .21938  3934. 

1 

We  conclude  our  discussion  of  G(y)  with  an  asymptotic  analysis 
of  I^(y)  as  defined  by  (18).   When  y  is  large  the  principal  contribution 
to  G(y)  is  from  I«(y)  and  is  of  the  form  (11).   When  y  is  small,  I„(y) 

is  small,  but  when  y  -  1,  I«(y)  is  of  the  same  order  as  I  (y) . 

2 
Consider  y  -  1  and  a  large  first.   Rewrite 

1  -2 

T  t    \  1     r   1  1   2  ,  /•    1-1    1  ,   2  ,  -a  log  y  dz    ,„t, 

12^^)  "  ^^P^ 2  ^^  ^^   J    exp{- jlog  z}z      ^  "   —    (31) 

a/2Tr       2a        0         2a 

2 
When  y  -  1  and  a     is  large  the  major  contribution  to  I„(y)  comes  from 

z  -  1  and  so  we  approximate 

i  ^  J_  .  1  _  (i_z)  +  (i_z)2  ~  2.  (32) 

z   z-z 

2 
If  we  replace  I/2  with  z  in  (31)  when  y  -  1  and  a      is  large,  we  see  that 

l2(y)  -  Ij^(y~''')  and  I  -  I  (y)  +  I  (y"""")  where  I^  is  given  by  (28)  and  (29). 

For  y  <<  1  we  write 

„    "  -2 

lo^y)  ^  ^^P''^ 2^°^   y^  /  exp{-z jios   2^^       ^       (^^^ 

a/27       2o        1         2a  z 


15 


2-2   2  2 

When  a   is  large  a   log  z  <<  z  for  1  £  z  £  exp{a  }  and  in  this  region 

1     2 
we  may  expand  exp{-  — j   log  z};  the  contribution  from  the  region 

2         2a  2 

z  >  exp{o  }  is   small    for  large  a     and  so  we  ignore  it.   Explicitly 

In(y)  =  -^  exp{-  -^log^y}  I     i^(-  -^)   -i—  r(a;l)     (34) 
o/27       2a       j=0  ^*    2a    9a  -■ 


where  for  negative  a,  r(a;l)  is  the  incomplete  gamma  function 


'^^•^'^  -  lohr  ^o'"'"  ^  ^''^ 


For  a  large  and  negative  a  standard  steepest  descent  calculation  yields 

vV     -e  -a         2  "^/^ 
^^"'^^  "ITO^^  ^   f2 +  e  +|]  (36) 

with  9  =  -(1  +  a  H — )  so  that  the  leading  term  in  an  asjnnptotic 
expansion  of  I„(y),  y  <<  1  is 

1         12         -2       -1       ' 

12^7)  -  ~  exp{ 2  l°g  y>[r(l  -  a   log  y)  ] 

a/2TT       2a 

,2^2      1/2  -  a"^log  y 

^2    log  y      ^2    log  y 

(37) 
2   2     -1/2 

By  differentiating  (37)  with  respect  to  ot,  higher  order  terms  may  be 
computed  if  desired. 
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2. 2   Summary  of  Approximations  to  G(y) 

2 
The  approximations  to  G(y)  for  a     large  computed  in  2.1  are: 

(1)  For  y  «  1  and  -(k+1)  £  ■^°^/   £  -k, 


2 

k      i       ° 

r^/  N  ~   V   (-y)     rl-2  2,  ,   1 ,   1  -   2  ,  ,  -2,     ,, 

G(y)  "  Z.  — TT  ^^Pi^yJ  ^  -f  ■*" e^Pi- olog  y)Y(o  log  y;l) 

j=0   ^*  a/27       2a 

k                 2                2 
,  (-1)   rl    .1,2,  ka.,^,l^,ka-,         1         , 
+  ^,       [-^expi — 2"1°S  (y^   )}erfc( ^log(ye   )) — ] 

2a  /Za  a/lTjlk+a      log  y] 

k+1  2  ? 

+  {^li..p(-i5loS^ye"=*"°   )}erfc(^los(ye "'«>''  ) 

2a  /2a 


]  +  io(y)  (38) 


a/2TT'[k+l+a      log  y] 

-2 

where  y(a     log  y;l)  is  an  incomplete  gamma  function  whose  singularities 

-2      -1 
at  integer  -k  and  -(k+1)  are  cancelled  by  the  singularities  of  [k  +  a  log  y] 

-2     -1 
and  [k  +  1  +  a  log  y]  respectively,  and  I^(y)  is  given  by  (37).   An 

asymptotic  evaluation  of  I„(y)  when  y  is  small  (cf.  formula  [37])  shows 

it  to  be  small  by  comparison  with  I^ (y) . 

(2)  For  y  -  1  and  y  >  1,  with  I, (y)  given  by  (28), 


G(y)  ^  I,(y)  +  I/y  ^--J--^xp{-  ^log^y}  I  y^li^-o'hog^Y 


2  2 

1        ,1   2,  r          J.    ,   1   ^            a,,l      £/li  -la,, 

-  -rexp{:ra    }  [y  erfc( log  ye     )   + -erfc( log  ye      )]  (39) 

a/2  y     a/2 

T      9    °°     T      1       2j+l 

+  1  -  exp{-  -^log^}  I     -i— (1°^) 
2a      j=0  r(|+j)  a/2 

(3)  For  y  =  1. 

G(l)  -  1  -  exp{|a^}erfc(— )  +^^—[1  +  Ei(-l)  -  y]  (40) 

where  y  =    .5772156649  and  Ei(-l)  =  .219383934, 
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(4)  For  y  -  1  and  y  <  1,  (39)  with  y  replaced  by  y   and  I^ (y)  given 
by  (29). 

(5)  For  y  >>  1, 

1         12-2            12-2 
G(y)  -  exp{ 2l°g  y^r(a  log  y)exp{ 2  ^    (^  log  v)  ^ 

a-Zln  la  la 

1-2  -2  2  -2 

X  [1 ^  i|j'(a  log  y)(l  -  a  ^   {a     log  y))] 

2a 

2 
Figure  1  displays  the  graph  of  G(y)  for  a     =3.0  and  y  =  0  and 

an  approximation  to  it  using  the  above  approximation  formulae.   The  solid 

line  is  drawn  through  points  computed  to  five  digits  accuarcy  by  numerical 

integration  of  (1) . 


(41) 
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lliJ, 
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N 
3.   Inversion  of  [G(y) ] 

The  density  h(K)  of  the  sum  of  N  lognormal  random  variables  is 

hm   =  ^     I     exp{Ky}[G(y)]%.  (42) 

The  major  contribution  to  this  integral  comes  from  the  region  Ky  =  0(1). 
Since  G(y)  has  a  different  form  for  each  of  the  regions  y  <<  1,  y  -  1, 

and  y  >>  1,  the  functional  form  of  h(K)  is  different  for  differing  orders 

2 
of  K.   For  each  of  the  following  cases  N  and  a  are  assumed  large. 

(1)   K  >  NM  .   For  this  case  we  approximate  G(y)  by  (38)  and  write 

~'^'^1        1         1    2 
G(y)  -  e    a(y)  +  ^!^-  exp{ y^og  y>b(y)  (43) 

a/l-n  2a 

with  a(0)  =  1,   The  number  of  terms  in  a(y)  is  determined  by  the  magnitude 

N 
of  K.   Approximating  [G(y)]   by  the  first  two  terms  of  the  binomial 

expansion  of  G(y)  expressed  as  (43), 

M^y 

[G(y)]^  -  exp{-NM^y}  [a^(y)  +— exp{-  ^logS}b(y)a^"^(y)  ] , 

av2TT       2a 

the  integral  (42)  is  approximated  by  a  sum  of  two  integrals.   If  a(y)  =   1 

the  first,  that  with  integrand  exp{(K  -  NM  )y}a  (y) ,  vanishes;  when 
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a(y)  -  1  for  y  -  0  it  is  exponentially  small.   Consequently  we 
approximate 

h(K)  ^  2^  .   -^—  /   exp{(K  -  [N-l]M^)y  -  -^log  y}b(y)a"   (y)dy 

exp{-  ^log^(K  -  [N-1]M  )} 
-  ^        2^  -i-  (44) 


0^  K  -  [N-1]M^ 

X  b({K  -  [N-llM^r-'-)a^'-'-({K  -  [N-llM^}"""-) 

Corrections  to  (44)  can  be  computed  in  a  straightforward  fashion. 
The  approximation  (44)  to  h(K)  shows  that  for  large  K  >  NM  , 

h(K)  is  lognonnal-like;  i.e.  the  leading  term  of  an  asymptotic  expansion 

2 
of  it  is  composed  of  a  three  parameter  (y=0,  a    ,  (N-l)M  )  lognormal 

density  with  argument  K  -  (N-l)M^,  times  a  correction. 

(2)   K  -  NM  ,   This  case  requires  specification  of  the  scales  of  N 

and  K  -  different  scales  give  different  answers.   In  particular,  for 

2  5  2 

N  =  0(exp{2o  })  and  K  =  0(exp{-^  }) ,  a(y)  may  be  approximated  by 

12  N 

exp{-yM  + -jVy  }.   Upon  expanding  [G(y)]  written  as  (43)  and 

integrating  term  by  term  we  obtain 

.In 

2    2   N         N     2 
f(K)  -  (2Tra  )    b'\K)exp{-  ^log  K} 

o 


N-1   „  -\  (K  -  (N-j)M^)^ 

+  M(-)(2ttV(N-j))    exp{ 2(jj_.)v } 

j=0  ^ 


(45) 


X  {(2Tra  )   exp{-  — ^log  ^^^qy^^ }  b(|j^  _  (n:j)mJ^^  ^ 
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(3)   1  <  K  <  (N-l)M^ .   This  case  may  be  treated  in  the  same  way  as 
case  (1)  with  (N-l)M  -  K  replacing  K  -  (N-l)M  : 

h(K)  ^  -^^[(N-l)M^  -  K]~-^exp{-  ^log^((N-l)M^  -  K)  } 

(46) 
b([(N-l)M^  -  K]"^  a^"^([(N-l)M^  -  K]"^ 


(4)   K  <<  1.   Use  the  expansion  (11)  of  G(y)  in  its  asymptotic  region: 


h(K)  -  Yj^  •  — —  /  exp{Ky 2l°8  y 2  ^  ^^     ^°^   y)r(a  log  y)dy 


(-^)^exp{-  -^log^K}^  (47) 

0^/27       2a^      ^ 


„  2.  ;:r---(N-l)^N,   2,    ^.    ,   N   ,2,   log  K,  , 
X  N   (a/2TT)     'r  (-a  log  K)exp{ r  i(j  ( 9")  ^ 

2a        a 
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